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We study the effects of quantum vortex fluctuations in two-dimensional superconductors using a dual 
theory of vortices, and investigate the relevance to underdoped cuprates where the superconductor- 
insulator transition (SIT) is possibly driven by quantum vortex proliferation. We find that a broad 
enough phase fluctuation regime may exist for experimental observation of the quantum vortex 
fluctuations near SIT in underdoped cuprates. We propose that this scenario can be tested via 
pair-tunneling experiments which measure the characteristic resonances in the zero-temperature 
pair-field susceptibility in the vortex-proliferated insulating phase. 



PACS numbers: 74.40. +k, 47.32. Cc, 74.50. +r, 74.76, 



I. INTRODUCTION 

The emergence of the pseudogap phase in cuprate su- 
perconducting materials may be attributed to strong or- 
der parameter phase fluctuations. Based on the small 
superfluid phase stiffness (SPS) at T = 0, together with 
its empirical scaling with the transition temperature (T c ), 
a proposal was put forward that the normal state pseu- 
dogap phase in the underdoped regime is a superconduc- 
tor whose phase coherence is destroyed by thermal phase 
fluctuations but with a robust gap-like feature due to 
a strong pairing amplitudeu a. Some experimental evi- 
dence was provided for strong thermal fluctuations of un- 
bound vortices over a widej-range of temperatures above 
T c in the pseudogap phasea. This suggests that the su- 
perconducting transition is of Kosterlitz-Thouless type 
with a broad phase fluctuation regime. 

It is then natural to explore the effects of strong 
T = quantum phase fluctuations, since the strengths 
of both thermal and quantum phase fluctuations are 
correlated with the small magnitude of SPS. For in- 
stance, quantum phase fluctuations reduce the Debye- 
Waller factor (e^) where <f> is the order parameter phase, 
and have significant consequences for the c-axis optical 
conductivityH, thejjenormalization of SPS, and the pair- 
field susceptibility^. A more dramatic effect of quan- 
tum phase fluctuations may be the quenching of phase- 
coherence by vortex pair proliferation in the underdoped 
regime, whichJeads to a superconductor-insulator tran- 
sition (SIT)Lra. Some experimental findings are con- 
sistent with- the existence of a quantum critical point 
near SIT&E2I, controlled by the charge carrier doping. 
However, direct evidence for vortex pair proliferation is 
yet to be found, although an experiment on current- 
voltage (J— V) characteristic in Bi2212 shows indications 
of a large density of quantum-, vortex pairs well be- 
low T c and far away from SITO. When searching for 
definitive experimental tests, it is important to construct 
and study reliable phenomenological theories which in- 
corporate vortex fluctuations. The critical properties 
of vortex-proliferated SIT have been studied using a 



framework of the dual transformation. In this frame 
work, dual vortex fields are conveniently introduced, 
the new order parameter of vortex proliferationQ'cS 
and the transition is described by the 2+1-dimensional 
(2+1D) Ginzburg-Landau theory of the quantum vor- 
tex order parameter, although the actual critical proper- 
ties are determined by the existence of long.-range in- 
teractions, periodic potential, and disordeiO. In this 
paper, we directly apply the mean-field dual formula- 
tion to two-dimensional (2D) superconductors and study 
its phenomenological relevance to cuprates in both su- 
perconducting and insulating states. We find that the 
phase fluctuation regime may be broad enough for exper- 
imental access. In vortex-proliferated insulating states, 
the characteristic form of the pair-correlation function is 
( e i*(*,t) e <*(o 1 o)} K e i °V / l*l 3 -l x l 2 /^ ) w ith a > 0. We pro- 
pose that this novel form of the pair-correlation function 
can be verified by pair-tunneling experiments. 



II. DUAL THEORY OF QUANTUM VORTICES 

In this paper, we denote the space-time 3- vectors with 
x and k and the spatial vectors with x and k. We use 
Greek indices running from to 2 and Roman indices 
i,j from 1 to 2. We also adopt a metric tensor g^ v — 
_ 2(5^°) to evaluate contractions of indices. For 
notational convenience, we set U = 1. We begin with a 
T = BCS model of a 2D superconductor coupled to the 
electromagnetic field. 
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where c CT is a fermion field, D„ = 
ant derivative, is the electromagnetic vector poten- 
tial, and m is the effective fermion mass. Here we do not 
specify the pairing symmetry, since it is not significant 
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within the accuracy of our discussion. Sf-A^] is an elec- 
tromagnetic gauge field action which takes the following 
form in 2D: 
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where we chose the Coulomb gauge V • A = and d is 
the thickness of the film. We assume that the order pa- 
rameter amplitude fluctuations are negligible [| A(x, t)\ — 
A], and focus on the phase degree of freedom <f> — 
— i ln[A(x, t)/A]. In order to decouple the (j> field from 
the order parameter amplitude, we perform a singular 
gauge transformation ^ CT (x,i) = c CT (x, £) e _l< ^ x '*)/ 2 , with 
iptr the field operators for the transformed quasiparticleO. 
We eventually arrive at the following effective theory of 
the phase and the electromagnetic fields after inte- 
grating out the fermion degrees of freedomO: 



(1) 



-(2A + V0) 2 ] +d t </>n a } +S[A^} , 



where p s is the SPS defined as p s — n s /4m, n s is the 
superconducting fermion density, and c s is an analog of 
the phonon velocity in a superfluid. In a 2D BCS su- 
perconductor, c s is related to the Fermi velocity vp hy 
c s = v F /y/2. We are interested in deriving the effec- 
tive theory of the phase fields which describes vortices. 
Therefore, we separate 9 M into d^tfi = 8^9 + A^ where 
9 is the spin-wave type Goldstone fluctuation and A^ 
gives topologically non-trivial phase gradients generated 
by vortices. Then we integrate out both 9 and A^ fields 
to obtain the final effective action of A u '- 



S[A,}^lY,^o(k)\A \ 2 ~IC T (k)\A\' 



(2) 



with a constraint that V-A = 0. Here /Co = ftk /(k 
u 2 + 8we 2 p s \k\) and K T = Ps{~uj 2 /c 2 + k 2 )/(A~ 2 + k 2 - 
w 2 /c 2 ) where A is the penetration depth. Over the time 
scale of our interest, JCt ~ p s k 2 /(A~ 2 + k 2 ). Below we 
will neglect the term d^cj) n s /2 in Eq. (|l|), which plays 
the role of a dual magnetic field but is insignificant near 
SIT since n s is renormalized and approaches zero in the 
insulating state. 

Now we assume a certain distribution of N point-like 
vortices. At a fixed time, the vortices in 2D superconduc- 
tors are pancake-like. In 2+1D, however, we can consider 
the space-time paths of pancake vortices as 3D lines of 
vortices. The current density (J7^) of the vortices are 
obtained from J» = e»» x d v A\. can be expressed 
in term of the space-time paths of vortices in 2+ ID as 
follows: 



1=1 ^ Ul 



ui) - x) , (3) 



where ui parameterizes the space-time path Xi of the l- 
th vortex. Therefore, Eq. (g) describes the interactions 
between infinitesimal vortex segments in 2+ ID (RefJla) if 
we re-express S'[,4 M ] in terms of J^. Here we assume only 
vortices of one flux quantum; the ones traveling backward 
(forward) in time are antivorticcs (vortices). 

The field theory of Eqs. (|J) and (^|) is inconvenient for 
description of particle-like vortices. Therefore, we first 
transform Eqs. (g) and (^) into particle dynamics. Later 
we will conveniently transform the particle dynamics into 
field dynamics of vortices. Before the transformation, we 
first separate the vortex-current interactions in Eq. (j^) 
into contact (short-range) and long-range interactions. 
The contact interaction can be approximately expressed 
in the form of a relativistic particle action as follows: 
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where we have deduced the rest mass from the contact 
interaction as following: 



m v vl $(°>(0) = (2^) 2 ^/C T (fc)/2k 2 



(5) 



The velocity v v is taken as v v = Cvf where C = 0(1) 
on physical groundsE3. In addition to Eq. (Q), there 
is contact repulsion between vortices of distinct labels. 
We will not estimate or discuss this repulsive interaction 
explicitly, except to mention that it stabilizes the non- 
zero expectation value of vortex fields in the insulating 
state which we will discuss later. In order to incorporate 
the long-range interactions into the particle dynamics, we 
first separate the dynamics of and A^ by introducing 
a Lagrange multiplier G^. This amounts to adding a 
term J d 2 x dt G^ vX d u A\ - J*) to the action in Eq. 
(^) which enforces the relation between and A^. Then 
we integrate out A^ fields to obtain the following effective 
interactions: 



k 



k 2 + A- 



a _ (k 2 c 2 + 8^e 2 p s |k|-^ 2 ) 



-2tt^ J dui G^ (Xi ) dX^/dui , 



(6) 



with a Coulomb gauge constraint V • G = 0. Thus, com- 
bining Eqs. (Q) and (||), we have obtained a theory of 
relativistic particles with a rest mass m v , coupled to the 
field Gfj_ which mediates the superfluid phase modes be- 
tween the vortices. 

We are now in a position to transform the particle 
dynamics into a field theory. Since we are interested 
in vortex-antivortex pair creation (annihilation), we now 
specialize into distribution of 2+ ID vortex loops by re- 
quiring that the paths X(u) are closed trajectories. Here 
we follow Ref£i3 and use the particle-field correspondence 
to re-express the total action S coa t + Stat in terms of a 
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relativistic complex scalar field $ coupled to the gauge 
field G M : 

S[$,$*,G] = J d 2 x dt \—(dn — 2niG fj)^* (d^ + 2niG^)<P 

-m$\$\ 2 ]+S[Gp] , (7) 

where G^ represents the local U(l) gauge symmetry of 
$ fields. Here S[G^] is the part quadratic in G M from 
Eq. (|) and we set w v = 1 for notational convenience. 
This is the well-known dual form of the theory of super- 
conductivity where the roles of magnetic (vortices) and 
electric (Cooper pairs) charges are interchanged^ 2 ! It is 
easy to show that the vortex current can be related to <f> 
as = i($9 M $* - $*<9 M $) + 47rG A1 |$| 2 , which au- 

tomatically satisfies the necessary conservation condition 

V = o. 

III. STRENGTH OF QUANTUM VORTEX 
FLUCTUATIONS IN CUPRATES 

Next, we wish to explore the feasibility of experimental 
observation by studying the width of the phase fluctua- 
tion regime where we expect to observe strong phase fluc- 
tuation effects. We assume a linear scaling of SPS with 
charge-carrier doping, and search for the magnitude of 
SPS (doping) at the most likely point of SIT in the under- 
doped cuprates. The action in Eq. (Q) can be viewed as 
the effective 3D Ginzburg-Landau functional of <I>. Then 
the SIT occurs upon ordering [U(l)-symmetry breaking] 
of the $ fields when to 2 < 0, so that ($) ^ In this 
case, ($(x) $*(?/)) w $o, where $q is a positive constant 
in the mean-field approximation, independent of \x — y\. 
The correlation function (<&(x)$*(y)) can be considered 
as the expectation value of the number of vortex paths 
that connect x and y. Therefore, ($(x)$* (y))= constant 
implies a non-zero and constant probability of arbitrarily 
long 2+1D vortex loops. The non-zero expectation value 
of $ leads to recovery of the U(l)-symmetry of the su- 
perconducting order parameter, and hence the Mcissncr 
effect is absent even in the presence of a non-zero bare 
SPS. Below we discuss the renormalization of m 2 to one 
loop expansion, assuming that the loop expansions are 
reliable. 

We use following parameters of optimally doped 
cuprates: A w 20 meV, A w 200 nm, the coherence 
length £ rs 2 nm, and the film thickness d rs 1.5 nm 
which roughly corresponds to a monolayer film. We as- 
sume a vortex core size of £ which provides a momentum 
(short-range) cutoff at A c = ir/£ and a frequency (short- 
time) cutoff at «i?A c . With this prescription, we avoid 
ultraviolet divergences and obtain 5^{0) = v F A 3 J{2n) 3 . 
We can then estimate m v and v v from Eq. (ra) in terms 
of the above parameters. The result is m v u^ « p s 7rln«; 
where k = AA C . We find that m v ~ 0.36 eV in opti- 
mally doped cuprates. For simplicity, we assume that the 
charge-carrier doping only affects magnitudes of the bare 




FIG. 1. (a)-(c) Diagrams of one-loop corrections to mj. 
The solid (dashed) lines denote $ (G M ) field propagators, (d) 
Solid (dashed) line is the one-loop-renormalized (bare) vortex 
rest mass (m v ) as a function of the bare superfluid phase 
stiffness (p s = n s /4m). m vo and p ao are respectively the 
vortex rest mass and the superfluid phase stiffness at optimal 
doping. 

SPS, and we keep the other parameters constant. From 
Eq. (Q) we can calculate one-loop corrections to m 2 as 
shown in Fig. |l| and find the point where the correction 
is of the same order of magnitudes as the bare value. 
In fact, the only correction comes from Fig. |l|(a); Fig. 
|l|(b) docs not contribute because the $ loop vanishes, 
and Fig. |l|(c) has been already included as the contact 
interaction. The correction 5m 2 /m 2 roughly behaves as 
— 8(i> v A c ) 3 /3pg7r 4 In 3 k, whose magnitude becomes large 
when the bare SPS is smaller. Eventually, we find that 
to 2 — (5m 2 — > (breakdown of Ginzburg-Landau theory) 
occurs when p s « 0.3A/hiK « 2 meV for cuprate films, 
which corresponds to about 10 % of that of the optimal 
doping [See Fig. [3(d)]. This implies that, using the em- 
pirical scaling between p s and T c , the phase fluctuation 
regime begins when T c is less than 10 K or so. In terms 
of charge-carrier (hole) doping concentration p, this cor- 
responds to p ~ 0.35p o where p is the optimal doping 
concentration. This is from the empiricaL-celation be- 
tween the critical temperature and dopingEZI, T c /T co ~ 
1 — 2.\(p/p Q — l) 2 , where T co is the optimal critical tem- 
perature and T c for p w 0.3p o . Therefore, there may 
be a broad enough doping range (0.3p o < p < 0.35p o ) for 
experimental studies on critical rVprtex fluctuations such 
as non-linear J—V characteristic!!!!. For instance, the un- 
derdoped YBCO used in the anomalous, proximity effect 
experiment performed by Decca et a/£3 fall within this 
parameter regime, which indicates strong phase fluctu- 
ations. The critical properties of the SIT are beyond 
the scope of this paper. Instead, we will simply assume a 
non-zero expectation value of and discuss definitive ex- 
perimental tests of vortex-proliferated insulating states. 
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IV. PAIR-FIELD SUSCEPTIBILITY 

In order to establish that underdoped cuprates are un- 
der strong quantum vortex fluctuations, a direct probe 
into the pair-fluctuations is .accessary. Here we propose 
a pair-tunneling cxperimentiiS to measure the pair-field 
susceptibility which contains information about pair- 
correlation functions. We first give a brief summary of 
the desired experimental setup. We consider a c-axis tun- 
nel junction of a thickness S between two cuprate sam- 
ples, where one of them is optimally doped with a c-axis 
penetration depth A c , and the other is an underdoped in- 
sulating film at T — with a thickness d. Since one of the 
electrodes is insulating, the usual Josephson current oscil- 
lating at a frequency of 2eV/Ti is absent. However, an ex- 
cess current will flow due to Josephson coupling of the su- 
perconducting pair-field of the superconducting electrode 
to the fluctuating pair-field of the insulator. Neglecting 
vortex fluctuations in the optimally doped electrode, the 
excess current can be related to the pair-field suscepti- 
bility of the insulator as I cx — (eE 2 / Sh 2 )lm D R {q,ui) 
where D R (x,t) = -i6»(t)([e^ xt \ e^' 00 )]) is the retarded 
pair-correlation function, Ej is the Josephson coupling 
energy of the junction, and S is the junction contact 
area. Here w is a frequency 2eV/h and q is a wave 
vector q(H) — 2eH(X c + d/2 + 5) /he which is deter- 
mined by a,-small magnetic field H applied parallel to 
the junctionlla. Therefore, the excess current can provide 
information about the spectrum of phase fluctuations. 
Below we obtain the form of the pair-correlation func- 
tion of vortex-proliferated insulators. 

First, we re-express the pair-correlation function as 
D(x,t) = — i(exp[i f*' „4 m (5;)g?5 m ]), ignoring contribu- 
tions from 8 fields since they are sub-leading in long 
length scales. Here we assume that 2eV < A to avoid the 
effect of order parameter amplitude fluctuations. Then 
we introduce a source function j^(x) = J**' dy^S^ (y—x) 
where we take y to be a straight line which connects (0, 0) 
and (x, t). Now we can re-express D as following: 



iD(x,t) 



f V<$>V®*VG, e l I ^^W+^t**** ,G M ] 

J M 

•£ fcM ^(fc)^(-fc)<^.(fc)^(-fe)>/2k 5 ' 



(8) 



where the second line is obtained by the method of func- 
tional integration (see Appendix |A|). In order to obtain 
the vortex-current correlation functions, we first define 
the polarization functions: 

U,(x - y)= {\<S>\ 2 8^{x -y)- V,(x)V;(y)/2). 

where V^x) = ^{x)d^*{x) - <^*(x)d^(x). We also 
define the longitudinal (Ho) and transverse (n^) polar- 
ization functions as the /z-component of the polarization 
functions where fj, is in the longitudinal (0) and trans- 
verse (T) direction respectively. Then the vortex-current 
correlation functions can be written as following 






2eV/E 

FIG. 2. Qualitative behavior of the excess current in an 
arbitrary unit. For comparison, the excess current due to the 
conventional finite-temperature fluctuations is shown in the 
inset. 
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[k 2 + 87r 2 /C T (£:)n o (£0]' 
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[k 2 + 87r 2 /Co(fc)nT(£0] ' 



(9) 



Using Eq. @ and assuming H^tc) « constant > 0, 
which holds in the vortex proliferated state, we obtain 
the asymptotic behavior of £)(x, t) in the large |x|, \t\: 



iD(x, t) cx exp 2ni 

in(L/c)n V^ 



(10) 



where L is the size of the sample. Upon Fourier trans- 
formation, we find that D{q,u) cx n /(cj 2 - E 2 ) 2 where 
E 2 = [27rhi(L/£)n T ] 2 +vlq 2 (H). Accordingly, the ex- 
cess current behaves as I cx (q,2eV) cx l/[(2eU) 2 — E 2 ] 3 , 
resembling the imaginary part of resonance peaks located 
near 2eV — ±2ir ln(L/£)Hr. The resonance peak heights 
are determined by the normal-state junction resistance 
and (| ( i ) | 2 ). The apparent logarithmic divergence of the 
exponent in Eq. (|l0|) and of the resonance energy is due 
to the fact that the energy of dual vortices (Cooper pairs) 
is logarithmically divergent as the system size, similar to 
vortices in superfluid helium. This weak divergence does 
not pose a serious problem in realistic samples which have 
finite sizes, however. This excess current qualitatively 
differs from that due to fluctuations at finite-temperature 
superconducting transitional!! where 7 ex cx uj/{lo 2 + T 2 ) 
(See Fig. ||). Therefore, it is possible to confirm the exis- 
tence of phase-fluctuation driven insulating states by de- 
tecting characteristic resonances in the pair-field suscep- 
tibility. The pair-correlation function oscillates at large 
\t\ due to the presence of vortex condensates. In the 
phase-coherent state of the vortex fields the trajecto- 
ries of Cooper pairs act as dual vortex paths in 2+1D, 
and the papjcorrelation function is the probability of a 
dual vortexN of length y/v^\t\ 2 - |x| 2 . Therefore D(x, t) 
is determined by the action of a vortex line which con- 
nects (0, 0) and (x, t). This explains the time-like length- 
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dependence of the exponent of -D(x, t). In fact, due to dis- 
sipative processes that we have not considered, the pair- 
correlation function decreases in magnitude at large \t\ in 
addition to the pure oscillation in \t\ shown in Eq. (10). 
Accordingly, the resonance peaks in the pair-correlation 
function will be broadened depending on the strength of 
the dissipation. 



S c = i d 6 x d*y(J A (x)Ji>(v))hi(x)h li (v)/2 . (Al) 



In the Coulomb gauge, the Fourier transform of Eq. ( Al) 
can be expressed as 

Sc =iJ2^( k )U-k){Mk)M-k)}/(2}< 2 ) . 

k . fx 



V. SUMMARY AND CONCLUSIONS 

We discussed the possibility of observing quantum vor- 
tex fluctuations in undcrdoped cuprate superconductors 
near SIT. Using the dual theory of vortices, we showed 
that cuprate superconductors are subject to strong vor- 
tex fluctuations so that it is possible to experimentally 
access the fluctuation regime near SIT. As a definitive 
test of the phase fluctuation scenario, we proposed an ex- 
periment to measure the pair-field susceptibility to probe 
the form of the pair-correlation function in the insulating 
regime. We expect that the pair-field susceptibility in 
vortex-proliferated insulating states is qualitatively dif- 
ferent from that observed in the normal-state fluctuation 
regime, and shows characteristic resonance peaks. The 
result can be generalized to any phase-fluctuation driven 
SIT of superconducting films. We anticipate that a more 
realistic dual theory of vortex fluctuations in supercon- 
ducting films can be obtained by including the effects of 
normal fluids or d-wave nodal quasiparticles. 
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APPENDIX A: PAIR-CORRELATION 
FUNCTIONS 

Here we give a brief derivation of Eq. (|J). From the 
first line of Eq. (||), the first term in the exponent can 
be rewritten as 



where h\ is defined through the relation j M = ie^ v ^d u h\. 
Assuming a Gaussian functional integral, we use the 
Baker-Haussdorf formula to obtain 

(J J d 3 xj x h x ^ _ gl S c _ 
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